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FIELD OF THE INVENTION 

The present invention relates to a method intended for real-time modelling, by 
neural networks, of hydrodynamic characteristics of multiphase flows in transient phase 
in pipes. 

The method finds applications notably for modelling of hydrocarbon flows in 
pipelines. 

BACKGROUND OF THE INVEOTION 

T™»pontag hydrocarbo.^ from pn>dac,ion sites u> tteating plants constitutes an 
important Unk in the pet«>leum chain. I. is a delicate Unk became of the complex 
interactions between the phases forming the transported effluents. Ihe basic objective 
for operators is to reach an optimum produOivity under U.e best safety conditions. They 
therefore have to control as best ti»y c«. the velocity ^rfflte t«nperatare so as to avad 
^necessary pressure drops, unwanted deposits and unsteady-state flows. The mea>od 
it^ is generally used consist, in modelling in fl.e best possible way the transportation of 
complex multiphase streams so as to provide a. all times an image of the flows m ib. 
various parts of the pnrfuction chain, taking mto account the precise constitution of the 

effluent, tiie flow rates, the pressures and the flow i^mes. 

There are currently various «>ftwa.» tools for simulating the transport of complex 

multiphase streams, allowing to design suitable production equipments at an early stage. 

Patents US-5.550.761. FR-2.756.044 (US.6.028.992) and FR-2,756,045 (US- 
5.960.187) filed by ti« applicant notably describe modellmg metirods and tools 
allowing to sunulate flte ttansport of complex multiphase streams on steady ot transient 
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flow and capable of taking into account instability phenomena that occur because of tite 
i„egular geometty of the formation crossed by the pipe o, of the topography thereof 
refeired to by specialists as « terrain slugging » or « severe slugging ». 

The simuUtion tools are as complex as the modelled phenomena. Precision and 
5 performances can only be obtained after, relativdy long modelling time, which is not 

really compatible with real-time management. 

AruAer approach aUowing, alone or m paraUel with the above modeUing methods, 
real-time m«««ement of the par«neters of a fluid circulation uses neural networks. 
I, can be teminded fl«t neural netwofts define a data processing mode simulating 
,0 the functiomng of biological neural systems. In such networks, an element carries out a 
relatively simple calculation «.ch as a weighted sum of the signals present at its inputs 
applied to a non-linear fenction. which determines flte state of the output fltereof. A 
large mmtber of such elements, interconnected in series and in paraUel, is used. Suiti*le 
selection of tite weighting fcctors aUows the network to carry out complex Actions. 
.5 Networks known as retropropagation networics for example use multiple layers of 
elements as defined above. Adaptation of such . networit » a precise task is carried out 
by "training- the network on a cettam mmiber of examples and by adjusting d» 
weighting &C«s for each element to ti-e suitable values. Input values are presented to 
.he networic. the output vatoe produced by tite network is analysed and fl» weighting 
20 factors are modified so as to minimis in tiie best possible way the diffe^nce betiveen 
fl« efi-ective value at the output and tiie value exp«ted in the example selected. After 
sufficient training. ti>e neti»ork is suited to respond to new input values ft.r which the 
output value is not known a priori and to produce a smtable output value. In its 
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pri.„=iple. a neural ne«™rk works according to a non-linear regression method which is 
aU fte more effective in relation to coiwentional methods. Two network types can be 
used, mainly MLP (MuM Layer Perceptron) networks, or Kohonen networks, well- 

known to specialists, 

5 Pat«« EP-1,176,481 filed by the appUcant describes a method intended for real- 
toe estimation of the flow regime, a. any point of a pipe whose structure is defined by 
. c«.»n mm*er of stn.«ure ^ physical parameters, of a multiphase fluid stream 
defined by several physical quantities and comprising liquid and gas phases. According 
to this method, the flow regime is modeUed by forming a non-linear neural networic 

10 with an input layer having as many inputs as there are stmoture parameters and phyrical 
quantities, an output layer with as many outputs as there are quantities necessary fi-r 

estimation of tite flow regime, and at l««t one intennediate layer, by forming a learning 
base wifl. predetermined tables comtecting various values obtained fcr the output data to 
the corresponding values of the input data, and by determining, by Herations, weighting 
,5 factors of flte activation iiincfion allowing tx> properly com«ct the values in tire input 

and output data tables. 

Output data of the neurons is preferably analysed so as to sort, among the values of 
the output data of the neural network, only the pertinent data to be taken into account 
for iterative detennination of the weighting coefficients of the activation fimction. 

20 Patent EP-1.217,474 also filed by the applicant describes a method allowing to 
construct a module (hydrodynamic or thermodynamic for example) that it is best suited 
to fixed operating conditions depending on the structure of the pipe and on a set of 
determined physical quantities (hydrodynamic or thermodynamic quantities for 
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^ fixed v»i«ion for *c pa«me.ers a=d ^ physical ,uan.Hies. 

The leaning base is adapted .o the imposed operating conditions ^ optinuzed neural 
oen^rks be^ «ljusKd « ti« imposed operating conditions are generated. In fl.e case, 
for example, where .he .^^le is to be integrated in a general multiphase flow 
,in.ul«ion modeU both hydrodynamic and tiaennodynamic. ti.e model is „«d to form 
the learmng base so as to select the set of physical ,uaMities ti^t is best suited to ti>e 
operation of tiae model, as well as ti.e variation ranges fixed for said parameters and s«d 
physical quantities, and tite optimized neural networks ti>a. best adjust to the learning 
base formed are generated. 
„ to the aforementioned prior m«hods. the flows are considered in a global way. 
^out distinction between flte various possible flow regimes of tite fluids in tite pipe : 
^ed flow, dispersed flow, intenmtten. flow, whose behaviours are different. This 
^ lead to modelling errors that are too great in relation to flte estimation -r-ality 
.e<,„ired for production momtoring. Furti«rm«e. ^ do not take account of tite 
.5 existence of simple models (for example analytic models) tr^lating in mathematical 
form characteristics of one or more flow r^es. 

SUMMARY OF THE INVENTION 

THe object of ti« method according to tiie invention is ti.e constiuction of a model 
for real-time simulation of flte hydrodynamic behaviour of a multiphase fluid flow in 
» transient phase in a pipe, considering fixed operating conditions concennng a certain 
number of determined structure parameters relative to a« pipe and a set of determined 
physical quantities with fixed variation ranges for said parameters and said physical 
quantities. Neural networks are used with inputs for sttucture parameters and physical 
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.^es. »d ou«„»s Where resuHs ™y ft-r of <^ ^V^o-V"-" 

behaviour are available, and a. lea. one intermedial layer. «.e neural ne^vorks being 
ae^r^ined i«™tive.y so a. «, adjus. .o .be vah.es of a .earring base wiU> predetennnred 
.bles connecting different vah.es obtained for ti. o«put dau .o the corresponding 
values of the input data. 

The method comprises : 
. constructing several neural networt. respectively dedic^edto differed 

re^mes, 

. a probabiUty neural network s«ted to evaWe a. all ti»es the 

, probabiUties for the flow in pipe to correspond respectively to the various flow 
regimes, and 

w the various neural networks weighted by said 
- combining the results provided by the vanous ncui 

probabilities. 

According to an implementation example, method comprises constructing a. 
„ lea. three neural networks respectively dedicated to the stratified flow regime, the 
dispersed flow regime and tire intermi.t«« flow regime, evaluating the probabiUties for 
^ fluid flow in the pipe to correspond respectively to tite ti»ee flow regimes and 
Unearly combining ti,e resuhs at ti. outputs of ti« ti^ee dedicated neural networl. hy 
weighting them by said probabilities. 

When the available database is sufficiently detailed to distinguish subregimes 
^titin a Single flow regime, a probability n«.ral networlc (RN^ suited to evaluate at 
any time the probabUities for ^ flow in the pipe «> corr^pond respectively to U« 
^ous flow subregimes distinguished in the various flow regimes is constructed and 
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*e re».« provided by v^ious .e»al o«wo*s .re combined b, weigbting 0^ by 
said probabilities. 

The estimation results obtained are aU the more accurate as : 
. a neurai mode, i, develops! by flow regime or ^.bregime. which aUows to take 
3 account of the parUcu.ari.ies of the physics contained in each one of the .aws 
represented, and 

. ^ cooti^rous and derivab.e (in the mathemaUca. sense) connection ...owing 
«^ition between «,e various .aws is cn=a.ed by . specific neura. or expert networic. 
Besides. ti» me^od retains the c.p.city of fl« tforementioned methods for 
,„ perfcrming re=d-time simulMion of the flows, and the results obtained benefit ftom the 
regularity of the estimation function obtained. 

BRIEF DESCRIPTION OF THE FIGURES 

Other features and advantages of the method according to the invention will be 
Clear from reading the description hereafter of a non-limitative embodiment example, 
15 with reference to the accompanying drawings wherein : 
. Figure 1 shows a model structure example, and 

. Figure2showsanexampleofthestructureofeachneuralnetworkofFi^el. 

DETAILED DESCRIPTION 

we consider . circu.ation of n^pbase fluids in a pipe with « .east a h,uid ph.se 
.„ .nd at .east a gas ph»e. and we try to construe, a mode, allowing, from a certain 
number of geometdca. and physic, input data reiative to the pipe »d of physic. d..a 
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of flow regime. We 4erefore us^ as me„tio««l above, for a given luanti^, S 
(Fig 1). various neural neworks specifically suited to various flow regimes in the 
pipe we construct for exan^le an expen network E» »odeUing specificaUy the 
«^ed flows, another one. modelling specfflcally the tatermittent flows and a 
one. modehing specifically the disperse flows. We also construa a neur^ 
mode, im^ specificaUy intended to evahrate a. time the probabiHt, p^ P. and 
p^. If S„ and are respecUvely the oun».t values of the three experts, wethen 

construct an evaluation function 5 such that : 

5 = Pstra Ss« + pDisp Soisp + Ps* Sm • 

Inputs and outputs of .h. various .«.r.l or e,p«. uhkW. forming the 

hydrodynamic model 

Whatever the neural model considered, the input data result from : 
. geometrical data : diameter of the pipe, roughness, inclination, etc., 
5 . data describing the characteristics ofthe fluid: density ofthe phases, viscosity 

the phases, etc., 

. a^chanu«rizingthennxture:gasflaction,gas«,uidsurftcetension.etc.. 
. linear or non-linear combinations of these inputs. 

. and also simplified models, continuous or not. containing informafion on the physics 
20 of the mixture. 

Each model produces for example at tire output the hydrodynamic behaviour of the 
effluents, and notably the flow regime. It evaluates and delivers at two mam outputs 
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velochy difference dV be«,een gas «.d liquid for exa,»ple. fte linear pressure drop 
aP/ac or the ft«*ion P » a [0 ;11) of flow of tte regime processed thereby. Other 
cmumties qualifying d>e flow regime can be calculated 6om these two oufl=uts. 

The outputs provided by the e^^erts are e,se«idly the velocity differences b«ween 
the phases u«ler the assumption of a certain flow regime (for example, the Stradfied 
expert delivers the esHmation of the velocity difference between the phases under the 
assumption of a stratified flow). 

m outputs provided by the probability network are the probability of belonging to 
each flow r^me processed by the «q>ert networks, knowing the taputs. 



Structure of the networks 

The various neural or expert networks dedicated to the different flow regimes are 
preferably Multi Layer Perceptron (MLP) type networks weU-known to the man skiUed 
in the art. generally estimating a hydrodynamic quantity. They comprise each (Fig.2) an 
„ input layer consisting of a certain number N, of neurons corresponding to the N, input 
dau of the complete physical model, an layer of a neuron for example 

corresponding to the param^er sough. (dV. aP/ac or p). and at leas, one intermediate 
layer, referred to as Mdden layer, whose number of n»n>ns N. is optimized. The 
^er of hidden layers and the number of neurons which constitute them are 
» ae.enmnedfh>mfl.e network leannng and validation resuKs. The network is totally 

connected, "nre non-linearity of tins network is obtained by either a sigmoid activ^ion 
function govenung the behaviour of the neun^s in the hidden layer, or by the identity 
fiinctionorthesoflmaxfunctionfortheoutputlayer. . 
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Tte neural networics comprise an input layer, one or two hidden layers and an 
output layer. The activation fonctions of the various neurons, well-known to the man 
skUled in the art. are either the sigmoid function (for the hidden layers), or the identity 
function or the softmax function (for the output layers). 
5 Learning 

The weigto of each network o, «per. are de«rmined « 4e end of a leamtog 
«age ; during .his sUge, .he are suppUed wi* a se. of dau forming .heir 

,e,™ng base, and *e config«ation and fte weigh., of Ae ne««* are optinuzed by 
n^nimizing erro. observed for all U,e s«np.es of fl.e base. be«,ee„ Ae o„.pu. da« 
.„ resuHing from ne«vork catadaUon and 4e daU «cpe«ed a. *e o««». given by 4e 
base. The erro. ean be fl« absolu« errors be^een *e inpu. and ou.p«. quantities or 
.he relative errors, according .o *e perfbmmce desired for ti.e ne^ork. Tire 
generalization capacities of ti,e ne.work are ti.=n «s.ed ftom i« capaci.y .o pn>periy 
ealculate tiie two outputs fcr inpuB to. are unknown fl^reto. 
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The databases used are of different naftjres : 
. for estinution of the velocity difference dV or of ti.e pressure drop, each base 
contains pairs of inpu./output values, each output value being ti.e desired value of 
^ estin^ted quantity in the case of the flow reginte processed by the dedic«ed 

network, 

,„ . for probability estin.ation.ti. desired output is a vec«>r of magmtude equal .oti» 
number of flow regimes considered (in tite example of Fi^e 1. ti>e vector is 
dimension 3); tins vector contains <N^ - D ««, values, and a value equal to I. 
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^ch ccn^ponds .o to pn*abi.hy for to fluid flow regime i« to pipe » 
correspond «> to flow regime processed by to dedicated neural ne«,ork. 
in to above e^pie. we have considers, toee different flow regimes : straSfied. 
i^er^i^ent and dispersed. «s is in no way iimirative. In cases where more detailed 
aata is available, allowing to mato distincdons within a single flow regime, for example 
. separate wavy stratified .on, smooth suatifled within to stra^fied fl- regime, 
specific experts modeUirtg each one of these subre^ - Prefer*>y created. 

Results 

By implementing such modeUin,, we obtain a continuous and infinite,, derivable 
^ient hydrodynamic model which carries out real-time calculation of to mam 

.ydrodynamic quantities characteri^ to flow. The pr^iry ^ 
.„ows.ocrea.eanov«-lhydrody.amic,aw.omtodifferen.flowlawsmodelledby 

^ various dedicated neural models. The transition between two flow laws is more or 

less steep (more or 



less great derivative) depending on to precision given to to 



, probabiUtyestimaUon. but it is continuous. wHcheUminates possible un^^^ 
^ ntodel results related to to existence of discontinuiries. The overall model is «..ed 
either use htdependentl, of any otor mo*Ue or integration in a complete model. 



